Draft version February 2, 2008 

Preprint typeset using lAT^^X style cmulateapj 



THE MACHO PROJECT LARGE MAGELLANIC CLOUD VARIABLE STAR INVENTORY. XL 
FREQUENCY ANALYSIS OF THE FUNDAMENTAL MODE RR LYRAE STARS 

Alcock, C.^'^'"*, Alves, D. R.'', Becker, A.**, Bennett, D.^''', Cook, K.H.^'"*, Drake, A.,^'^°, 
Freeman, K.-^^, Geha, M.-^, Griest, K}'^, Kovacs, G.^, Lehner, M.^, Marshall, S.^, 
MiNNiTi, D.^", Nelson, C.^'^^, Peterson, B.", Popowski, P.^**, Pratt, M.^'', Quinn, P.^*^, 

RODGERS, A7, StUBBS, C*'^^, SUTHERLAND, W/^, VANDEHEI, T }^ AND WELCH, D.L/* 

Draft version February 2, 2008 
ABSTRACT 

We have frequency analyzed 6391 variables classified earlier as fundamental mode RR Lyrae (RRO) 
stars in the MACHO database on the Large Magellanic Cloud (LMC). The overwhelming majority (i.e., 
96%) of these variables have been proved to be indeed RRO stars, whereas the remaining ones have fallen 
in one of the following categories: single- and double-mode Cepheids, binaries, first overtone and double- 
mode RR Lyrae stars and non-classified variables. Special attention has been paid to the properties of 
the amplitude- and phase- modulated RRO stars (the Blazhko stars). We found altogether 731 Blazhko 
variables showing either a doublet or an equidistant triplet pattern at the main pulsation component 
in their frequency spectra. This sample overwhelmingly exceeds the number of Blazhko stars known 
in all other systems combined. The incidence rate of the Blazhko variables among the fundamental 
mode RR Lyrae stars in the LMC is 11.9%, which is three times higher than their rate among the first 
overtone RR Lyrae stars. No difference is found in the average brightness between the single-mode and 
Blazhko variables. However, the latter ones show a somewhat lower degree of skewness in their average 
light curves and a concomitant lower total amplitude in their modulation-free light curves. From the 
frequency spectra we found that variables with larger modulation amplitudes at the higher frequency 
side of the main pulsation component are three times more frequent than the ones showing the opposite 
amplitude pattern. A search made for a modulation component with the Blazhko period in the average 
brightness of the individual variables showed the existence of such a modulation with an overall amplitude 
of sa 0.006 mag. On the other hand, a similar search for quadruple modulation patterns around the main 
pulsation component have failed to clearly detect such components at the ~ 0.004 mag level. This means 
that the amplitudes of the quadruple components (if they exist) should be, on the average, at least ten 
times smaller than those of the triplet components. This finding and the existence of Blazhko variables 
with highly asymmetric modulation amplitudes not only question the validity of the magnetic oblique 
rotator model, but also puts stringent constraints on models based on mode coupling theories. 

Subject headings: globular clusters: general — stars: horizontal-branch — stars: oscillations — stars: 
variables: other (RR Lyrae) 



1. INTRODUCTION 

The mystery of the physical cause of amplitude- 
modulation in certain RR Lyrae stars, popularly known 
as the "Blazhko Effect" remains despite almost a century 



of study. In a previous paper in this series (Alcock et 
al. 2000, hereafter AOO) we discussed the results of the 
frequency analysis of a sample of 1350 Large Magellanic 
Cloud (LMC) variables, formerly classified as first over- 
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tone RR Lyrac (RRl) stars. In addition to the discov- 
ery of more than one hundred double-mode variables, we 
found 52 stars which displayed long-period amplitude and 
phase modulations. Until very recently (Olech, Kaluzny & 
Thompson 2001, and references therein; Cseresnjes 2001; 
Moskahk & Poretti 2003; Soszynski et al. 2003) it was 
established as a well-observed phenomenon only among 
fundamental- mode variables. The study of AOO gave the 
first reliable statistics on the incidence rate of this behavior 
among RRl stars. A rate of 4% was found which is sur- 
prisingly low compared with the most often referenced rate 
of 20-30% for the fundamental mode (RRO) stars (Szeidl 
1988). The main purpose of the present study is to ver- 
ify whether this large difference does indeed exist between 
these two groups of stars. The large and homogeneous 
sample of the macho database stars used in this study 
yield much more significant results for the global statisti- 
cal properties of this population than any other existing 
investigation which were typically based on very limited 
and inhomogcneous samples of objects from various stel- 
lar systems. (Exceptions to this characterization are the 
surveys of Soszynski et al. 2003 on the 7612 RR Lyrae 
stars of the LMC from the OGLE-ii database, of Moskalik 
& Poretti 2003 for the 215 variables of the Galactic Bulge 
from the OGLE-i database and of Cseresnjes 2001 for the 
3700 variables of the Fornax Dwarf galaxy.) 

Apart from the task of identification of the variable type 
for RR Lyrae in the LMC MACHO database, we also study 
the global properties of the Blazhko stars and compare 
tlicin with their singly-periodic counterparts. A key dis- 
criminant between rival explanations for the Blazhko be- 
havior is found in the predicted frequency and amplitude 
pattern distributions of the amplitude spectra of these 
stars. The existence (or non-existence) of triplet or quin- 
tuplet structures and the degree of asymmetry of the cor- 
responding modulation amplitudes strongly restricts the 
range of acceptable models. The present large dataset 
enables \is to put significant constraints on the possible 
ranges of the various modulation components. A brief 
summary of the results obtained by this analysis of a frac- 
tion of the present dataset has already been published by 
Welch et al. (2002). 

2. THE DATA, THE METHOD OF ANALYSIS AND 
VARIABLE CLASSIFICATION 

The dataset analyzed in this paper comprises 6391 light 
curves selected from the MACHO database as variables sat- 
isfying certain conditions on brightness, color and period, 
corresponding to the expected properties of fundamental 
mode RR Lyrae stars in the LMC^. Due to overlapping 
fields, 364 light curves turned out to represent variables 
already appearing in the database with other identifica- 
tion numbers. Since these double (or multiple) identifica- 
tions decrease the number of stars in each variable type by 
only a few percent, and therefore, does not infiuence the 
incidence rates in a significant way, all of our subsequent 

^The selection of the RR Lyrae for this paper was made many years before the completion of the MACHO Project. In the time since the 
submission of this paper, we have become aware that there axe no RR Lyrae variables in the eastern half of field 5 and an underabundance in 
a small number of other locations. This was apparently due to an access problem of the program used to phase the photometry. Our ability 
to test this hypothesis end-to-end has been hampered by the unavailabiUty of scripts present on Mount Stromlo Observatory computers which 
were lost during the devastating fire of Junuary 18, 2003. We have since extracted a more complete list of RR Lyra« from the full eight-year 
database and we will make this list available to interested investigators. It is our strong expectation that the frequency of Blaahko behavior 
reported in this paper will not be significantly altered by a similar analysis of the additional candidates in the eight-year list. 

■^http; / /www.macho.mcmaster.ca 



results and statistics refer to the full sample of 6391 iden- 
tifications. The list of stars with multiple identifications 
is presented at the end of this section. 

This paper analyzes RR Lyrae stars identified in 30 
LMC MACHO Project fields (1, 2, 3, 5, 6, 7, 9, 10, 11, 12, 13, 
14, 15, 17, 18, 19, 22, 23, 24, 47, 53, 55, 57, 76, 77, 78, 79, 
80, 81, 82). These fields include 54% of the observed stars 
in all fields because they cover the densely populated bar of 
the LMC (see the home page^ of the macho Project). The 
data utilized span a period of 2700 d (7.4 yr), which is 1 yr 
longer than the timebase of the first overtone RR Lyrae 
stars analyzed in AOO. The magnitudes have been trans- 
formed to the Johnson V system according to the method 
described by Alcock et al. (1999). As a result of the longer 
time baseline, the average number of observation epochs 
per star is also greater, often in the range 700 to 1000. 
The sampling rate and data distribution is the same as in 
the case of the RRl stars. In particular, there is only a 
1 d quasi-periodicity in the sampling but no trace of 1 yr 
aliasing. The accuracy of the individual data points is also 
the same as for the RRl data, but, because of the larger 
amplitudes, the overall signal-to-noise ratio (SNR=A/cr, 
A is the total (peak-to-peak) amplitude, a is the standard 
deviation of the residuals) is larger ((SNR) « 11, whereas 
for the RRl stars (SNR) « 2.G, although this figure refers 
to the data in the 'r' color, which further lowers the SNR 
- see AOO). 

Apart from spectrum averaging (sec Sect. 5), the 
method of analysis in this survey is essentially the same 
as that in AOO. We performed a standard frequency anal- 
ysis based on a Fourier approach and a 'pre- whitening' 
technique (e.g.. Deeming 1975, with the fast transforma- 
tion method recommended by Kurtz 1985). Before execut- 
ing the Fourier transformation, outlying points have been 
omitted by checking the distribution of the magnitudes 
and leaving out data points separated from a compact set 
which is thought to represent the lightcurve plus random 
noise. The relative number of points omittc;(l in this way 
varied between 1% and 4% and most often these omitted 
points were obtained in very poor seeing. The classifica- 
tion of the variables was performed in three major steps: 

• Frequency analysis with pre-whitening in the 
[1,6] d~^ band, visual inspection of the frequency 
spectra and folded lightcurves, and preliminary se- 
lection of variable types. 

• Reanalysis of the pre-selectcd multiperiodic variables 
by three successive pre-whitenings and a closer ex- 
amination of the frequency spectra. 

• Visual inspection of the lightcurves of all variables 
to check the above classifications. 

Except for the first pre-whitening cycle, when the main 
pulsation component is filtered out from the time series, 
we used a single-component least squares fit for whitening. 
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Table 1 

Variable types in the macho LMC database 



Type 


Short description 


Number 


% 




RRO-S 


Smgly-pcnodic b U RR Lyrac 


4882 


79.3 


0.5 


■p> T~) n TOT 1 

RKU— rsLl 


RRL) witn 1 close irequcncy component 


400 


0.5 


0.3 


n F> A TOT O 

RRU— BLi 


RRO with 2 close symmetric frequency components 


331 


5.4 


0.3 


RRO-i/M 


RRO with several close components 


20 


0.3 


0.1 


RRO-PC 


RRO with period change 


177 


2.9 


0.2 


RRO-D 


RRO & integer frequencies 


308 


5.0 


0.3 


RRO-MI 


RRO with some miscellany 


40 


0.6 


0.1 


RRl 


Singly-periodic FO RR Lyrae 


36 






RROl 


FU/FO RR Lyrac 


6 






NC 


Non-classific;d variables 


32 






SOCEP 


Singly-periodic overtone Cepheid 


60 






FU/FO 


FU/FO double-mode Cepheid 


2 






FO/SO 


FO/SO double-mode Cepheid 


58 






BI 


Eclipsing binary 


39 







Note. — The symbol <j{%) denotes the standard deviation of the population ratio, assuming 
Poisson distribution (see Equation (1)). In computing population ratios, only fundamental 
mode RR Lyrae stars have been considered. 



Because of the high nonlinearity, for the main pulsation 
component we utilized all harmonics up to the third one. 

It is important to emphasize that the classification is 
based on the observed frequency pattern in a wide fre- 
quency band. Therefore, the significance of a given pattern 
is determined through a comparison of a large number of 
other features which might make that pattern less signifi- 
cant than if we made the comparison in a smaller frequency 
band (the chance of detecting higher, noise-induced peaks 
increases with the increase of the frequency band). For 
example, if we choose a narrow band around the main 
pulsation component, and search for equidistant triplet in 
the spectrum, a considerable fraction of the Blazhko (BL) 
stars that were classified after the inspection of the wide 
frequency band as BL stars with doublet patterns (here- 
after BLl stars) will prove to have triplet patterns and will 
be classified as BL2 stars (see Sect. 4). Although some ar- 
guments can be given for the preference of a classification 
based on pattern search in a narrow frequency band, we 
follow a more conservative; approach, and our classification 
(if not stated otherwise) refer to the results obtained from 
the above wide frequency band. 

Variable types and their respective totals obtained in the 
course of the above analysis are listed in Table 1. The no- 
tation is the same as in AOO, except for the Blazhko stars, 
where we use similar symbols both for variables with ap- 
parently single modulation components (BLl stars) and 
those with two symmetrically-spaced modiilation com- 
ponents straddling the main pulsation component (BL2 
stars). In AOO we assigned the symbol "BL" only to the 
BL2 stars, in order to maintain a coherent nomenclature 
with the handful of variables which have been analyzed in 
the Galactic field. These few variables all display triplets 



in their frequency spectra. However the analysis of our 
larger sample strongly suggests that the two types of mod- 
ulation overlap and they are indicative of the same under- 
lying phenomenon which we denote by "BL" . 

As discussed earlier, our statistics have not been cor- 
rected for multiple identifications due to field overlap. This 
results in incidence rates which differ from the true rates 
by, at most, 0.1%. One reason that we did not make this 
correction is that, in critical cases, classification of the 
same variable star may depend on the particular time se- 
ries we analyzed. For example, the same star might appear 
in one field with many observed epochs and in the other 
with a much smaller number. If a modulation component 
has an amplitude near the noise level, it could potentially 
be classified as both a Blazhko and as a single-mode vari- 
able (an example is with variable star identified as both, 
MACHO ID 17.3194.3230 and 10.3194.499). A similar am- 
biguity might arise in the classification of BLl and BL2 
stars if the signal-to-noise ratio is low (e.g., macho ID 
6.6811.969 and 13.6811.4172). We foimd five variables 
which have been classified both as BLl- and BL2-types. 
Among the 400 BLl stars there are 13 duplicates with 
the same BLl classification. Similarly, we found 8 dupli- 
cate identifications among the 331 BL2 stars which were 
classified identically in both instances. These duplications 
suggest that our classification is repeatable and reliable. 

The notation in Table 1 incorporates the notation of 
Udalski ct al. (1999) for the radial pulsation mode, in 
which the symbols FU, FO and SO denote the hmdamen- 
tal, first and second overtone modes, respectively. The 
RRO— z/M variables are similar to the BL stars in that 
their amplitude spectra exhibit close secondary compo- 
nents. However, their frequency patterns are more com- 
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plicated and show no obvious sign of dominant, equidistant 
frequency components (around the main pulsation compo- 
nent). For PC variables, the structure remaining around 
the main pulsation component cannot be resolved, and the 
few successive pre-whitenings do not eliminate these fea- 
tures from the amplitude spectrum. Although the failure 
of the pre-whitening technique indicates a long-term varia- 
tion in the signal frequency (or even in the amplitude), we 
must note that it is also possible (at least in some cases) 
that we are witnessing Blazhko modulations with very long 
periods. This contention is supported by the fact that the 
present sample contains BL stars with modulation periods 
of several years. With this long BL periods, the modu- 
lation components are quite close to each other but are 
still well-resolved due to the long temporal baseline of the 
observations. Several BL stars additionally display a PC 
character. In such cases the variables have been classified 
as BL-type, despite the additional PC behavior. 

A significant fraction (5%) of the stars in this sample 
display frequency components at integer d^^. These are 
labelled RRO— D variables (such stars also were found in 
AOO). Clearly this behavior indicates some instrumental 
or reduction artefact. At present the source of this is un- 
known. 

The difference between the RRO— MI and NC variables 
is that in the former class (most probably) all variables 
are RRO stars, but they exhibit extra peaks in their fre- 
quency spectra, which are not easy to explain within any 
simple framework (e.g., multimode pulsation). In the case 
of NC stars, classification is not possible either because 
of the limited dataset, or because of the 'strange-looking' 
light curve (e.g., drift-like mean-light variation). While 
additional interesting physical behavior may be present in 
these stars, we expect (but have not established) that ex- 
trinsic effects are responsible for many such classifications. 

Many new LMC Cepheid variables have already been 
discovered, both by the MACHO and by the OGLE Projects. 
After cross-identification of the apparent Cepheid vari- 
ables found in our initial sample with both databases, we 
find that 67% of the FO/SO Ccpheids reported in this 
paper are new discoveries. 

Statistical errors (standard deviations cr(%)) of the in- 
cidence rates are calculated from the assumption that the 
population follows a Poisson distribution (i.e., = \/Wi, 
where Ni is the number of objects in the i-th population). 
This yields the following expression for the standard devi- 
ation of the population ratio Ni/N 

^^■/^ = ;^/W-¥) • 

BL stars constitute a significant fraction of the RRO pop- 
ulation in the LMC according to Table 1. Therefore, for 
the first time, BL stars can be used for reliable statisti- 
cal studies. These results will be presented in the next 
sections. 

Since the spatial distribution of variables is a separate 
issue from the ones to be discussed in the following sec- 
tions, here, in Figure 1, we show the field- by- field distri- 
bution of the BL and all RRO stars. It is seen that there 
are considerable differences in the number of RRO stars 
among the various fields, but this correlates with the stel- 
lar density of the fields (e.g., fields 53, 55 and 57 are far 



from the LMC bar). The ratio of BL to total number 
of RRO stars also shows some fluctuations, especially in 
the lower-density fields. Therefore, this fluctuation might 
have a statistical origin. For example, for a sample size of 
N=50, an observed incidence rate of 3% has a 2.4% stan- 
dard deviation. In the well-populated fields the relative 
number ratios fluctuate around 10%, suggesting that this 
ratio might be the same in all fields. When the full sample 
is considered, we get an average incidence rate of BL stars 
of (11.9 ± 0.4)% (the error corresponds to la standard 
deviation of the mean) . 




MACHO Field 



Fig. 1. — Field-by-field distribution of the number of all funda- 
mental mode RR Lyrae stars (lower panel, with the field number 
shown on the top of the columns) and the ratio of the Blazhko stars 
to the total number of the fundamental mode RR Lyrae stars (upper 
panel). 

We recall that the incidence rate of the RRl— BL stars 
in the LMC is 3.9% (see AOO). Although the frequency of 
the RRO— BL stars have proven to be lower than the com- 
monly cited value of 20-30%, they are still three times 
more common in the LMC than their RRl counterparts. 
It is possible that some RRl-BL stars might have been 
missed because of their low modulation amplitudes but it 
seems unlikely that so many BL stars were selectively lost 
due to the lower signal-to-noise ratio. 

Tabic 2 lists the identification numbers, classification 
and comments concerning all individual frequency spec- 
tra. Primary periods (those which appear with the high- 
est amplitude in the frequency spectra of the original 
data) are listed in Table 3. Fourier decompositions of 
all BL stars are given in Table 4. Denoting the funda- 
mental mode and Blazhko (or modulation) frequencies by 
fo and fm, respectively, a Fourier sum with frequencies 

— Jo ~ fm, V2 = /o, = /o + /m, = 2/o — /m, 

= 2/o, Vq = 2/o + fm., = 3/o, = 4/o, t'g = 5/o, 

vva = 6/o has been fitted to each light curve and out- 
liers were successively omitted until no data points devi- 
ated from the fitted curve by more than 3ct (i.e., three 
times the standard deviation of the residuals). The prin- 
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Table 2 

Star-by-star comments on the frequency analysis 



MACHO ID 


TYPE 


REMARK 


1.3441.1031 






1.3442.1051 






1.3442.1107 






1.3442.1243 


BL2 


? Long mod. period, Df=0.00063 (1/T=0.00037) 


1.3442.503 






1.3443.1313 






1.3443.1335 






1.3444.1187 


RRl 


& PC 


1.3444.56 


SOCEP 




1.3446.1877 







Note. — This table is available in its entirety only in the electronic version of this paper. 



Table 3 

Co-ordinates and main periods of the variables 



MACHO ID 


a(J2000.G) 


(5(J2000.0) 


PERIOD [d] 


1.3443.1335 


05:01:26.0 


-69:15:38 


0.4509916 


1.3444.1187 


05:01:55.4 


-69:11:31 


0.4239698 


1.3444.56 


05:01:49.9 


-69:12:46 


0.8057676 


1.3446.1877 


05:01:42.7 


-69:02:06 


0.5997655 


1.3446.1886 


05:01:35.1 


-69:02:06 


0.6061866 


1.3446.2004 


05:01:37.9 


-69:03:26 


0.4876528 


1.3449.1327 


05:01:42.3 


-68:52:06 


0.4839085 



Note. — This table is available in its entirety only in the electronic version of this paper. 



Table 4 

Fourier decompositions of the Blazhko variables 



MACHO ID 


/o 


fni 


N 


V 


a 














«3 


0,4 


a--, 




0,7 






a 10 


rl 


y'2 


■f-'> 








'r~ 








1.3446.2004 


2.050641 


0.011805 


762 


19.1238 


0.0922 










0.0234 


0.2982 


0.0461 


0.0246 


0.1305 


0.0557 


0.0808 


0.0391 


0.0240 


0.0041 


3.7147 


1.7243 


4.0731 


1.1233 


5.6287 


2.2739 


3.5659 


1.6217 


5.9261 


2.8276 


1.3570.930 


1.709248 


0.008804 


304 


19.0981 


0.0750 










0.0138 


0.3227 


0.0641 


0.0284 


0.1512 


0.0255 


0.1033 


0.0613 


0.0365 


0.0174 


5.5545 


6.0614 


0.1049 


2.3508 


1.9917 


2.1388 


4.3896 


0.4926 


2.6433 


5.2757 



Note. — This table is available in its entirety only in the electronic version of this paper. 
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Table 5 

Variables with multiple identification numbers 



MACHO ID 


PERIOD [d] 


Q!(J2000.0) 


^(J2000.0) 


18.3446.4944 


0.6061888 


05:01:35.0 


-69:02:06 


1.3446.1886 


0.6061866 


05:01:35.1 


-69:02:06 


1.3446.2004 


0.4876528 


05:01:37.9 


-69:03:26 


18.3446.4648 


0.4876523 


05:01:37.8 


-69:03:26 


18.3449.4465 


0.4939059 


05:01:21.3 


-68:51:30 


1.3449.1459 


0.4939031 


05:01:21.3 


-68:51:31 



Note. — This table is available in its entirety only in 
the electronic version of this paper. 



Table 6 

Average Fourier parameters of the singly-periodic RRO and Blazhko stars 



Parameter 


Average(RRO) 


Average (BL) 


^0 


19.340 ±0.008 


19.351 ±0.010 


Ai 


0.293 ±0.003 


0.295 ±0.003 




0.136 ±0.001 


0.123 ±0.001 


^3 


0.094 ± 0.001 


0.073 ±0.001 


Atot 


0.843 ±0.008 


0.755 ±0.007 


(t>21 


2.389 ±0.005 


2.365 ±0.008 


(t>31 


5.076 ±0.011 


4.922 ±0.014 


(t>4,l 


1.571 ±0.019 


1.301 ±0.025 



Note. — The symbol Atot denotes the to- 
tal (peak-to-peak) amplitude without modula- 
tion. Errors show the standard deviations of 
the means as estimated from the scatter of the 
observed values. A 4cr-cut is applied in the com- 
putation of the Aq averages. 
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ciple effect of these omissions is to remove photometry 
acquired in conditions of very poor seeing. Amplitudes 
and phases refer to the following type of decomposition: 
V{t) = V + aisin{2'K{t — t^jui + ipi) + .... Here t is given 
in HJD— 2400000.0 and the epoch to is set equal to zero 
for all variables. 

Modulation frequencies, as listed in Table 4, have been 
computed through the combination of the visual inspec- 
tion of the frequency spectra and the automatic search 
for the "strongest secondary component" as described in 
Sect. 4. The visual inspection has resulted the revision of 
the fm values obtained by the automatic search for 19 vari- 
ables. The failure of the automatic search in these cases 
was caused by various sources, mostly by large power resid- 
uals at the main pulsation component, but high noise level 
and aliasing were also matters of concern. 

Multiple identifications due to field overlaps are listed 
in Tabic 5. For full time series of all variables we refer 
to the MACHO Project home page. Note that the time se- 
ries available there are for the full span of MACHO Project 
observations, unlike the data analyzed in this paper. 

3. COMPARISON OF THE GLOBAL PROPERTIES OF THE 
RRO AND BL STARS 

The purpose of this section is to compare the over- 
all photometric parameters of singly-periodic and Blazhko 
variables. We include all BL stars, because a preliminary 
investigation did not find any differences between the two 
subgroups with respect to the relations to be discussed in 
this section. 




.4 .5 .6 .7 

Period [d] 



Fig. 2. — Comparison of the average magnitudes V for the singly- 
periodic RRO stars (open circles) and for the RRO-BL stars (filled 
circles). 

First, we examine the average magnitudes. These val- 
ues are computed from a Fourier fit described in the pre- 
vious section. We used lOth-order fits for both the BL 
and singly-periodic variables. For the BL stars the fits 
include the basic pulsation component and its harmonics 
only up to order 6 because of the addition of the modula- 
tion components at /o and 2/o. We found that this fitting 
procedure gives reliable values for the low-order modula- 
tion and Fourier components. The average magnitudes are 
plotted in Figure 2. In order to avoid unnecessary crowd- 



ing, we took a subsample of 1000 variables from MACHO 
fields 1, 10, 11, 12, 13, 14 to represent the RRO-S stars. 

First of all, it is seen that both groups of stars cover 
the same ranges of parameters in this diagram. Based 
solely on their absolute magnitudes, Blazhko variables are 
indistinguishable from singly-periodic RRO stars. Indeed, 
as shown in Table 6, the average of the mean magnitudes 
of both classes agree within their statistical errors. Both 
classes of stars display large scatter in V. This is due to 
various effects, including crowding, inhomogeneous red- 
dening within the LMC, differences in the absolute mag- 
nitudes, finite geometrical extension of the LMC and, of 
course, measurement/calibration errors. Despite the large 
scatter, some trend in the average magnitudes is observ- 
able. This trend is in qualitative agreement with the well- 
known relation between period and luminosity observed 
in other stellar clusters (e.g., Kovacs & Walker 2001), and 
also expected from stellar evolution theory, which predicts 
brighter horizontal branches at lower metallicities (i.e., at 
longer periods - see Castellani, Chieffi & Pulone 1991). 
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Fig. 3. — Comparison of the Fourier amplitudes of the /o, 2/o and 
3/o components of the singly-periodic RRO stars (open circles) and 
of the BL stars (filled circles). 

Next we turn to the comparison of the Fourier param- 
eters related to the main pulsation component /o and its 
harmonics. Denoting by Ai the amplitude and by the 
phase of the Fourier component at frequency ifo, and em- 
ploying the standard, epoch-independent phase difference 
0ii = 01 — Figures 3 and 4 show the variations of these 
quantities as functions of the period. It is seen that with 
increasing order, the Fourier amplitudes of the BL stars 
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become significantly lower than those of the RRO-S stars 
(at the same period). A similar, but much weaker effect 
also exists among the phases as can be inferred from the 
computed average values and from their errors (see Ta- 
ble 6). This result implies that the modulation-free light 
curves of the BL stars show less skewness and have lower 
total amplitudes than their singly-periodic counterparts 
(see Table 6). The less significant difference between the 
phases suggests that the (Poj^ai) ^[Fe/H] relation de- 
rived by Jurcsik & Kovacs (1996) may have a wider va- 
lidity - when applied to the average light curve - than 
was originally suggested by these authors on the basis of 
the temporal Fourier decompositions of a few BL stars. 
(Lack of similarity between the temporal lightcurves of 
the BL stars and those of the single-mode RRO stars has 
recently been studied by Jurcsik, Benko & Szeidl 2002.) It 
is also worth noting that the definite correlation between 
the period and the phases implies that the metallicity of 
the majority of the RRO stars in the LMC is confined to a 
relatively narrow interval, as follows from the Pq ^ Pi/ Pq 
diagram for the double-mode variables (Popielski, Dziem- 
bowski & Cassisi 2000; Kovacs 2001), and also from direct 
spectroscopic measurements (Bragaglia et al. 2001). 
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Fig. 4. — Comparison of the relative Fourier phases of the 2/o, 3/o 
and 4/o components of the singly-periodic RRO stars (open circles) 
and of the RRO-BL stars (filled circles). 

4. MODULATION FREQUENCY AND AMPLITUDE 
DISTRIBUTIONS 

In this section we will discuss both the frequency with 
which modulation timescales appear in the BL popula- 
tion and the incidence with which the stronger modula- 



tion component is at higher or lower frequency than the 
primary pulsation peak. For this purpose, we define the 
sign of the modulation frequency /,„ (|/m| = |/o - /±|) to 
be positive if A+ > A- and negative if the reverse is true 
(here and A- are, respectively, the amplitudes of the 
modulation components on the higher and lower frequency 
sides of the main pulsation component). With the above 
definition we can include all BL stars and study if there 
is a preference for a dominant frequency/amplitude pat- 
tern. Figure 5 shows the resulting distribution function. 
It is clearly seen that there is a strong preference for fre- 
quency patterns where the modulation amplitude on the 
higher frequency side is larger than its lower frequency 
counterpart. We find that 74% of all BL stars have pos- 
itive /,„. Similar ratios are observed when the two types 
of BL stars are studied separately: /m > for 79% of 
the BLl stars and 68% of the BL2 stars. Blazhko vari- 
ables with very long modulation periods and those with 
periods as short as 50 days occur with about the same 
probability. For shorter modulation periods the number 
of BL stars sharply decreases and reaches very small in- 
cidence rates for modulation periods shorter than about 
25 d. Nevertheless, one can still find a few stars with very 
short Blazhko periods: 9.5360.804 and 6.6212.1121 have 
fm = 0.050 and 0.122 d~^, respectively, with the latter 
value being the only one in the whole sample exceeding 
0.1 d^^. We also note that similarly short modulation 
period was observed in AH Cam by Smith ct al. (1994). 
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Fig. 5. — Distribution of the modulation frequencies for the 731 
BL stars. See text for the sign convention of fm- 

It is also a matter of interest whether the small devia- 
tions from a strictly equidistant frequency spacing in the 
case of the BL2 stars are statistically significant. We recall 
that in the course of a similar analysis of the RRl stars 
it was shown that the observed deviations are not signifi- 
cant (see AOO). Following the same notation as in AOO, we 
compute the quantity Sf = f+ + f- — 2/o to characterize 
the degree of deviation from equidistant spacing. We show 
the distribution function obtained from the analysis of the 
BL2 stars in Figure 6. Wc select the triplets from the fre- 
quency spectra trough the following automated procedure: 
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• Compute successively pre-whitened amplitude spec- 
tra in the ±0.1 d^^ neighborhood of the main pulsa- 
tion component. Five pre- whitening cycles are com- 
puted and the highest peak obtained after each cycle 
is stored. 

• The "strongest secondary component" / is cho- 
sen from the stored peak frequencies by requiring 
\f — fo\ > ci/T, where ci = 1.5, T is the total times- 
pan of the respective dataset. This condition pro- 
duces a distinction between PC and long-period BL 
stars. 

• A third peak is selected from frequencies appear- 
ing on the side opposite to the main component 
where the "strongest secondary component" was 
found. The variable is classified as a BL2-type if 
\Sf\ < C2/T, where C2 = 3.0. This condition is used 
to select triplets with "approximately" equidistant 
frequency spacings. 

In nine cases out of the 331 BL2 variables, this automatic 
selection method could not find a triplet, usually because 
the large PC component coexisting with the triplet did not 
allow the process to reach the level of the Blazhko com- 
ponents within the five pre-whitening cycles. A high noise 
level was also a factor in some cases. Therefore, Figure 6 
shows the results for the remaining 322 stars. (It is noted 
that five variables from this sample yielded considerably 
smaller /„ values than the ones listed in Table 4. Again, 
this difference is accounted for by the large PC component 
in these stars. In order to maintain a close compatibility 
with the subsequent statistical tests, we keep these dif- 
ferent modulation frequencies in the investigation of the 
distribution of the non-equidistant spacings.) 

We see that, with few exceptions, all frequency devia- 
tions are within ±0.0002 d~^. This is about one half of 
the characteristic linewidth (i.e., l/T w 0.00037 d^^). In 
order to assess if this degree of unequal spacing still can 
be explained as a result of noise-induced frequency-shift, 
we perform the following test. Artificial time-series are 
generated with equidistant triplets given by the Fourier 
decompositions (see Sect. 2) of the 331 BL2 stars. Then, 
a Gaussian white noise is added to these data with the 
standard deviation of the fit of the original data. In the 
next step these test data are processed in the same way as 
the real observations. Finally, we use the peak frequencies 
obtained from the analysis to compute Sf, and plot the 
corresponding histogram. The result is shown in Figure 7. 
Comparison with Figure 6 reveals that the two distribu- 
tions are very similar to each other, although the test data 
show a smaller dispersion and no outliers above \Sf\ > 
0.0002 d~^. The derived standard deviations are consis- 
tent with this conclusion. After applying a 3<t criterion for 

6.3 X 10-5 d-i 

By performing the above 
test for 10 different realizations of the noise for each time- 
(4.59 ± 0.26) X 10-5 d-\ where the 

val- 
ues obtained with the various realizations. We conclude 
that the difference between the observed and test values 
is statistically significant. Therefore, the observed small 
deviations from equidistant spacing cannot be accounted 
for by purely noise-induced frequency-shifts. 



the omission of the outliers, we find a^j 
and aff = 4.9 x 10"^ 

10 difff 
series, we get a^y* 

error represents the standard deviation of the ten tr^^''* 



In checking the notes (see Table 2) related to the nine 
variables mentioned before for not passing the automatic 
selection criterion and those 18 stars which did not pass 
the 3(7 criterion with respect to the symmetric spacing, we 
find the following. Except for two cases, the lightcurves 
of the variables have some peculiarity - often of a PC na- 
ture - but we also find remarks on high noise level and 
small second modulation component. In one case there is 
no comment and in a single case the observed frequency is 
declared 'nice'. 
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Fig. 6. — Distribution of the differences in the triplet spacings for 
the BL2 stars. 
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Fig. 7. — Distribution of the differences in the triplet spacings 
obtained from test data as described in the text. 

Our conclusion from all these tests is that the cause of 
non-equidistant frequency spacing has various origins, in- 
cluding cases when neither the noise level, nor peculiar 
additional features are suspected to be responsible for the 
deviation. The cause of the non-equidistant spacing should 
be studied in detail on a star-by-star basis by considering 
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various effects in addition to the noise-induced frequency- 
shift as discussed above. 

In some models a resonance between the fundamental 
and Blazhko periods Pbl = l//m could be a matter of 
interest. Therefore, we compute the fractional part of 
Pb l/Pq and plot its distribution function. As can be ver- 
ified from Figure 8, there are no resonances between the 
two periods and the distribution is essentially flat for all 
fractional values. The fluctuations are of statistical origin, 
because the relative standard deviation in an average bin 
is 23% (assuming a Poisson distribution) . 
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Fig. 8. — Distribution of the fractional part of Pbl /Pa for the 731 
BL stars. 
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Fig. 9. — Distribution of the ratio of the maximum modulation 
amplitude to the main pulsation amplitude for the 731 BL stars. 

We now investigate the distribution of the modulation 
amplitudes. There are two important properties to ad- 
dress: (a) range of the relative (pulsation vs. modulation) 
strength of the modulation amplitudes, and (b) degree of 
asymmetry of the left- (A-) and right- (A+) modulation 



components. We employ the Fourier decompositions both 
for the BLl and for the BL2 stars to get the observed val- 
ues of Ai , A- and A^ . Since noise always introduces some 
non-zero excess power at all frequencies, a distortion of the 
distributions at low modulation amplitudes is expected. 
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We plot the distribution of the relative modulation 
strength in Figure 9. In order to minimize the effect of 
noise mentioned above, we take the largest modulation 
component in each case. The distribution sharply peaks 
at « 0.15 and decays with a long tail after « 0.2. It is 
important to observe that there is a wide range of possible 
modulation strengths from very weak to very strong with 
the largest having A±/Ai > 0.4. Considering extreme 
modulation strengths, we call attention to the specific 
cases of variables 2.5870.4967 with MAx(yl±/^i)= 0.55 
and 9.4997.754, with MAX{A±/Ai)= 0.10. The modula- 
tion components are clearly detected for both stars. 

The asymmetry in the modulation amplitudes is a signif- 
icant problem for both presently available models (Shiba- 
hashi 2000, Nowakowski & Dziembowski 2001). To quan- 
titize the degree of asymmetry, we introduce the following 
parameter 



A+-A^ 



A, 



(2) 



We recall that the amplitudes are computed from 
Fourier fits which include both the A- and A^ compo- 
nents. Q is found to cover almost all values, partly as a 
result of noise and partly due to the presence of 'hidden' 
BL2 variables among the BLl stars. In fact, if we perform 
an automated search for symmetric frequency patterns in 
the way mentioned at the beginning of this section, we 
find that 160 out of the 400 BLl stars could be classified 
as BL2-type. 

Figure 10 displays the distribution function of all BL 
stars. As we have already seen in Figure 5, there is a 
strong dominance of variables with A^ > A^. The dis- 
tribution is very broad but the most probable cases are 
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those with Q ^ 0.1 0.6. Of course, when the above dis- 
tribution is plotted only for the BLl variables, we find a 
deficit at small \Q\ values and the distribution is double- 
peaked at larger /smaller Q values. The BL2 variables are 
mostly confined to the Q = (—0.2,0.4) interval, but some 
extreme asymmetric cases can also be found. For example, 
variable 77.7793.1432 shows a large positive asymmetry 
with Q = 0.823. 

5. A SEARCH FOR ADDITIONAL MODULATION 
COMPONENTS BY SPECTRUM-AVERAGING 

Our standard method of analysis described in Sect. 2 

did not find any definite signature of two interesting fea- 
tures (both of these are expected within the frameworks of 
cither the oblique rotator or of the resonant mode coupling 
models). The first feature we expect is due to the long- 
term periodic variation of the average brightness. The 
presence of this component is expected from non-linear 
mode-coupling. The other possible feature is an equidis- 
tant quintuplet structure around /q. This property is pre- 
dicted by the oblique rotator model and is also expected 
when mode-coupling occurs. Since such features are not 
visible in the observed (individual) frequency spectra, we 
employ a method of spectrum averaging in order to am- 
plify such small signals (if they are present). 

Our procedure is the following. The huge differences 
in the lengths of the Blazhko periods require us to nor- 
malize the spectra in both the amplitude and frequency 
domains before averaging. Therefore, after pre-whitening 
by the 'trivial' components (i.e., fo - fm, fo, fo + fm, 
etc.) and taking the Fourier transform of the residuals, 
each spectrum is divided by its variance for more effec- 
tive noise filtering of the different quality spectra. (We 
always use amplitude spectra - both in the computation 
of the individual and of the average spectra.) Then, all 
peaks are identified and shifted by the amount necessary 
to position the modulation components in the same place 
in each case. Then, this unified frequency band is divided 
into equidistant bins and the maximum amplitude values 
of the peaks are found for each bin. These values are used 
in the binned-and-normalized version of original frequency 
spectrum. Finally, all these spectra are added together and 
the average is plotted for visual inspection. To demon- 
strate the effectiveness of the method, Figure 11 reveals 
the result obtained by the analysis of a small sample of 
BL2 variables. 

We see that the signal-to-noise ratio of the average spec- 
trum has dramatically improved compared to the major- 
ity of the individual spectra. The pre-whitening by the 
fo component (and its harmonics), results in a dip at 
Af / fmod = 0. As we shall sec in the; following subsections, 
dips related to pre-whitening will show up more obviously 
when using a large set of variables. Occasionally, a rather 
large separation is observable between the peaks in the 
individual normalized spectra (e.g., variable 5.4288.4400). 
This effect is due to the long modulation periods of those 
stars, which results in a greater stretch in the frequency 
domain. Therefore, because of the finite temporal baseline 
of the observations, we observe fewer peaks in the very nar- 
row frequency band corresponding to the long modulation 
period. PC behavior is also observable in some cases (e.g., 
MACHO ID 5.4285.302). The variable 5.4891.855 represents 
one of the handful cases where a quadruple structure is ob- 



served, but with unequal frequency spacings. 
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Fig. 11. — An example of spectrum-averaging. Individual normal- 
ized spectra are shown in each panel, except for the upper right one, 
where the resulting average spectrum is displayed. Eax;h spectrum 
has been computed after pre-whitening by fo and its harmonics. 
Only the peak amplitudes are displayed. The symbol A/ denotes 
/ ~ /O) where / is the test frequency. All variables are classified as 
BL2. 

5.1. Average light level 

In order to reach the highest possible signal-to-noise ra- 
tio, we employ all BL stars in the computation of an aver- 
age spectrum. Individual frequency spectra are computed 
in the [0,0.2] d~^ band. By using 200 bins in all cases, 
each spectrum is transformed in the frequency domain of 
< A///„, < 3, where A/ denotes the test frequency 
/. The result is somewhat sensitive to the choice of num- 
ber of bins. If the bins are too wide (compared to the 
linewidth), then the hidden signal will be less effectively 
recovered because of the neighboring smaller peaks (which 
come from the noise) will lower the amount of power col- 
lected in the bin of interest. On the other hand, bins which 
are too narrow will result in a division of signal power into 
several adjacent bins due to the noise- induced frequency- 
shift of the peak associated with the signal. With the 
above choice of bin number we cover approximately the 
characteristic linewidth at the moderately short Blazhko 
period of ~ 50 days. This bin width is also appropriate for 
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possible frequency shifts due to noise (see Figures 6 and 
7). In order to avoid too large a deformation in the plot 
due to the large dip at zero frequency (each time series 
has zero average), the final spectrum is shifted downward 
by an appropriate amount, then stretched in the vertical 
axis to cover the [0, 1] interval by the [min,max] range of 
the spectrum. 
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Fig. 12. — Average frequency spectrum of all BL stars in the 
proximity of the expected long-term modulation with the Blazhko 
frequency fmod- The symbol A/ denotes the test frequency /. Am- 
plitudes and relative frequencies are normalized as described in the 
text. 
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Fig. 13. — As in Figure 12 but a sinusoidal tracer with an 
amplitude of 0.010 mag and frequency 2f^od added to the signals. 

Figure 12 shows the resulting amplitude spectrum. We 
observe a definite peak at l^f / fmod = 1-0. Some ear- 
lier studies on Galactic Blazhko variables also indicated 
the presence of this component (Borkowski 1980; Kovacs 
1995; Nagy 1998). However, the data used in those studies 
suffered from very strong aliasing by the main periods of 
the variables caused by the observational technique em- 
ployed when those variables were observed. As a result, 
the possibility of power leakage from the /q ± /,„ compo- 
nents to frequencies differing by /o from these modulation 
components is expected. 

Since the observed amplitude is influenced by the size of 
bins and noise level, the simplest way to get an estimate 
on the amplitude and on the significance of the peak is to 
add a sinusoidal tracer of a given frequency to the signal. 
This has been done, and the result shown in Figure 13 is 
obtained. By comparing the two peaks we estimate the 
overall amplitude of the observed modulation in the aver- 
age light level of about 0.006 mag. 



A further characteristic feature observable in the spec- 
trum is the slow rise of the signal power from l^f / fmod ~ 
1.0 to zero. This mild trend is most probably due to long- 
term instrumental or photometric calibration effects. 

5.2. Quintuplet structure 

Side peaks at / = /o ± 2fm have been searched for in 
the same manner as the fm component. Here the number 
of bins used was 400 because of the larger bandwidth re- 
quired by this analysis. The average spectrum of all the 
BL stars is shown in Figure 14. In addition to some excess 
power very close to the prewhitened /o component and to 
/o±/m, we may also suspect a peak at /o — 2/m, very close 
to the noise level. The excess power near /o (and probably 
also at fozt fm) is associated with the observed PC char- 
acter of some variables. On average, this is a small effect, 
but in a few individual cases this may be comparable to 
the Blazhko effect (e.g. 6.6570.771 and 78.5615.905). 
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Fig. 14. — Average frequency spectrum of all BL stars in the 
neighborhood of /q. The main pulsation (/o) and modulation {/_ 
and /+) components have been prewhitened before Fourier trans- 
formation. The symbol A/ denotes f — fo, where / is the test fre- 
quency. Amplitudes and relative frequencies are normalized using a 
procedure described in the text. 




- ' " 



-3-2-1012 

A///™d 



Fig. 15. — As in Figure 14 but with a sinusoidal tracer of 

amplitude of 0.007 mag added to the signals at /^f / fmod = i2.5. 

In order to estimate an upper limit for possible quadru- 
plet components, we injected tracer signals into the indi- 
vidual time series in the same fashion as in our analysis 
for the fm component. The effect of an 0.007 mag ampli- 
tude tracer is shown in Figure 15. We see that this type 
of component would have been easily detectable from the 
present database. A test with an even lower amplitude of 
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0.005 mag shows that this would be also detectable, al- 
though with a lower significance. Therefore, we can safely 
conclude that if equidistant quadruplet structures exist in 
the BL stars, they would have amplitudes (on average) less 
than 0.004 mag. 

6. CONCLUSIONS 

We have studied the global statistical properties 
of amplitude- and phase-varying fundamental mode 
RR Lyrae (BL) stars by performing a frequency analy- 
sis of 6391 variables from the MACHO LMC database. We 
found 731 BL stars which is by far the largest popula- 
tion known in a single stellar system. Because of the large 
sample and the long time baseline of the MACHO Project 
observations, it was possible to obtain statistically signif- 
icant results, describing certain overall properties of the 
BL phenomenon. The main results are the following: 

• Fundamental mode (RRO) BL stars constitute 11.9% 

of the presently analyzed LMC RRO population 
(6158 stars). This incidence rate is three times larger 
than that of their first overtone (RRl) counterparts 
(see Alcock et al. 2000). 

• A continuous transition seems to exist between vari- 
ables showing an equidistant triplet around (and in- 
cluding) their main pulsation component in their 
frequency spectra and those displaying only a close 
doublet (BL2- and BLl-types, respectively). This 
suggests that strongly asymmetric modulation am- 
plitudes are not rare among the BL2 stars and their 
BLl coimterparts represent only extreme cases when 
the missing component has such a low amplitude 
that it becomes buried by the observational noise. 

• 74% of the BL stars have larger modulation ampli- 
tudes on the higher frequency side of the main pul- 
sation component. 

• Statistical tests indicate that some fraction of the 
BL2 stars may have non-equidistant frequency spac- 
ing. There is also a small siibgroup of 20 variables 
(RRO— i^M class) which display several close peaks at 
the main pulsation component but without apparent 
regularity in their frequency spacing. 

• There is a nearly uniform distribution of Blazhko 
modulation periods (Pbl) from « 50 to several hun- 
dred or even thousand days. In comparison, vari- 
ables with Pbl < 30 d are rare. 

• Modulation amplitude and A^) components 
around the main pulsation component (Ai) in the 
frequency spectra are usually in the range 0.1 < 
A±/Ai < 0.3, but there are also more extreme val- 
ues. 

• There is a variation with the overall amplitude of 
« 0.006 mag with a period of Pbl in the average 

brightness of the BL stars. 

• No clear trace of a quintuplet structure is found 
around the main pulsation component. The overall 
upper limit of the amplitude of such a component is 
estimated to be less than 0.004 mag. 



• No difference is found in the long-term average 
brightness betwcK;n the singly-periodic and BL RRO 
stars. However there are small, but systematic, 
differences in the Fourier decompositions of the 
modulation-free light curves of the BL stars and 
those of the singly-periodic stars. The differences 
yield average BL light curves which have smaller 
amplitudes and lower degree of skewness than their 
singly-periodic counterparts. 

• The PC variables (stars with unresolved power at 
the main pulsation component) show an incidence 
rate of 2.9%, which is much lower than the rate of 
10.6% of the RRl-PC stars. 

Most of the above conclusions are very similar to those 
given by Moskalik & Poretti (2003) for the Galactic bulge 
RR Lyrae stars. The main difference is that they find a 
much higher incidence rate for RRO-BL stars. Their rate 
is 23%, which is in the often-quoted range for the Blazhko 
stars (Szeidl 1988), but a factor of two higher than what 
we obtained for our LMC sample. Since the 215 RRO stars 
analyzed by Moskalik & Poretti (2003) can be considered 
as a statistically significant sample, we tend to agree with 
them that metallicity is the most probable agent for caus- 
ing the difference in the incidence rates. 

Very recently Soszyhski et al. (2003) published a sim- 
ilar analyisis on 7612 RR Lyrae stars in the LMC from 
the OGLE-n database. Although a detailed comparison of 
our study with theirs is beyond the scope of the present 
paper, we note the following. Their incidence rates for the 
BL stars arc 15% and 6% for the RRO and RRl stars, re- 
spectively. Both of these rates are somewhat larger than 
ours. It is unclear at this moment what is the reason of 
this slight discrepancy. It may come either from the higher 
signal-to-noise ratio obtained from the use of Differential 
Image Analysis, or from the higher density of variables in 
the LMC bar of the OGLE sample, or from the details in 
the search for the modulation components (for example, 
a search in a narrower frequency band around the main 
pulsation component may lead to higher detection rates - 
see Sect. 2). Further comparative studies are needed to 
answer these questions. 

Unfortunately, neither the oblique magnetic rotator 
(Shibahashi 2000) nor the rotating resonant pulsator mod- 
els (Van Hoolst, Dziembowski & Kawaler 1998; Dziem- 
bowski & Cassisi 1999; Nowakowski & Dziembowski 2001) 
are capable of explaining the observed properties of the 
BL population. One major problem for these models is 
the production of strongly asymmetric modulation com- 
ponents. No hint of this behavior is found in the oblique 
magnetic rotator model, whereas the resonant one is able 
to explain only mild asymmetry. Although no current ob- 
servational effort is known to be attempting to detect mag- 
netic fields in RR Lyrae stars, the last result in this field 
by Chadid (2001) does not support the presence of strong 
magnetic field in RR Lyrae itself. The situation is also 
complicated concerning the direct detection of nonradial 
pulsation through surface velocity field (Kolenberg 2002). 

We conchide that the underlying reason for Blazhko be- 
havior is still not understood. The large sample reported 
here will allow investigators to select subsamples with com- 
mon features for further study and will hopefully lead to a 
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theoretical understanding of this rather obvious observa- 
tional phenomenon. Precise spectroscopic observations to 
disentangle the large radial pulsation velocity field and the 
small nonradial one, together with a definite measurement 
of an upper limit on the magnetic field, are among the 
most important observational tasks for the near future. 
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